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ABSTRACT 

Extending the theory we derived recently for IID209458b to different cases of strongly irradi- 
ated gaseous exoplanets, we have calculated the consistent evolution of the new transiting planet, 
OGLE-TR-56b, for its recently revised mass determination. The theory is shown to successfully re- 
produce the observed radius, for the proper age of the system. We also examine the dissipation of 
kinetic energy at the planet's internal adiabat due to atmospheric winds, and place constraints on 
the efficiency of this process. We show that a fraction ~ 0.1 - 0.5% of the incident flux transformed 
into thermal energy deposited at the adiabatic level can accommodate the observed radii of both 
OGLE-TR-56b and HD209458b. The present theory yields quantitative predictions on the evolution 
of the emergent spectrum and fundamental properties of hot-jupiters. The predictions for radius, lumi- 
nosity, temperature as a function of the planet's mass and orbital distance can be used as benchmarks 
for future detections of transit planets. 

Subject headings: binaries: eclipsing - planetary systems - stars; individual (OGLE-TR-56, HD 209458) 



1. INTRODUCTION 

\ Since the discovery of 51 Peg b (Mayor & Queloz 1995) , 
over a hundred gaseous planets have been discovered in 
orbit around G, K, and M stars. These discoveries have 
opened a new domain of research in astronomy, at the 
crossroad between stellar and planetary physics. Many 
ground-based and space-based missions are planned for 
the coming years to search for more of these objects, and 
the next decade will probably see the first direct detec- 
tion of an extrasolar planet. These projects require a 
robust theoretical background to accurately understand 
the properties of these objects and to provide reliable 
guidance for observations. The vast majority of these 
planets have a projected mass, m x sini, of about one 
Jupiter mass (Mj„p ~ lO'^^Mo) but are found to or- 
bit much closer to their parent star than Jupiter, with 
a pile-up around a ^ 0.04 AU. As such, these planets 

' experience strong irradiation from their parent star. Ob- 
servations of objects that transit their parent star pro- 
vide information crucial to understanding these so-called 
hot-jupiters. The transit light curve or photometry, sup- 
plemented by radial velocity measurements of the star, 
yields the absolute mass and average radius of the planet. 
Together with an estimate of the age of the parent star 
from stellar evolution models, the radius and mass pro- 
vide stringent constraints on the evolution of the internal 
and atmospheric properties of strongly irradiated plan- 
ets. The correct description of these properties has been 
questioned recently by the detection of the transit planet 
HD209458b (Brown et al. 2001). No consistent model 
can adequately reproduce the observationally determined 
mass and radius of the planet {Rp = 1.42^0 ^3 Rj^p) 
(Cody & Sasselov 2002) without stretching these deter- 
minations beyond their error bars (Guillot & Showman 
2002, Baraffe et al. 2003, Burrows et al. 2003). This led 
to the suggestion that an important part of fundamental 
physics might be missing in the description of irradiated 
planets and has motivated several authors to propose 



possible shortcomings of the theory. In the present let- 
ter, the theory developed in Baraffe et al. (2003) for 
non-irradiated and irradiated planets is applied to vari- 
ous cases of hot-jupiters, with special focus on the case of 
the second observed transit, OGLE-TR-56b. The theory 
adequately reproduces the observed radius of the planet, 
for the proper age of the system and the more recent mass 
determination (Torres et al. 2003). The comparison of 
the present predictions for the evolution of hot-jupiters 
with forthcoming transit detections will allow the unam- 
biguous verification or rejection of the present theory. 

2. EFFECT OF IRRADIATION 

For a planet in close orbit around its parent star, the 
strong incident stellar radiation substantially modifies 
not only the emergent spectrum but also the evolution 
of the irradiated planet, as shown initially by Guillot et 
al. (1996). These effects are examined below. 

2.1. Atmosphere 

Under strong irradiation conditions, the solution of 
the transfer equation, which determines the atmospheric 
thermal structure of the planet, must include the incom- 
ing incident stellar flux in the source function. Such cal- 
culations were first conducted for hot-jupiters by Sea- 
ger & Sasselov (1998), using a fimited set of opacities, 
and have subsequently been improved by Barman, AUard 
& Hauschildt (2001) and Sudarsky, Burrows & Hubeny 
(2003). Figure n compares the thermal structures of ir- 
radiated and non-irradiated planets with surface grav- 
ity log 5 = 3.5 (cgs units) and intrinsic temperatures 
Toff — 500 K and Tcs = 100 K, respectively, orbiting 
a G2 main sequence star at 0.023 AU or 0.046 AU. The 
atmosphere models were produced using the same cloud- 
free opacity setup described by AUard et al. (2001) and 
Barman, AUard & Hauschildt (2001). As examined be- 
low, this Tcff sequence corresponds to an age sequence of 
-3x10^ yr to -9x10^ yr for OGLE-TR-56b. The irra- 
diated thermal profile is shown to be strongly modified 
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with respect to the non-irradiated one. As stressed in 
Baraffe et al. (2003), the differences between irradiated 
and non-irradiated temperature profiles clearly demon- 
strate the necessity to do consistent calculations that 
combine the irradiated atmosphere profile and interior 
profile. Matching arbitrarily the inner profile to an at- 
mospheric profile defined by the equilibrium temperature 
of the planet yields severely flawed results. 

Note that the present calculations assume that all 
the incident flux is concentrated on the day-side of 
the planet. Such an assumption maximizes the effect 
of irradiation and might overestimate the flux deposit. 
The present nearly isothermal irradiated profiles, how- 
ever, are consistent with the 2^m observed spectrum of 
HD209458b (Richardson, Deming & Seager 2003). Fu- 
ture calculations will include the latitude dependence of 
the incoming radiation but it is certainly interesting to 
verify whether the present calculations, with minimal as- 
sumptions, adequately reproduce the radius of new ob- 
served transiting planets. 

2.2. Evolution 

Consistent evolution calculations of irradiated plan- 
ets have been conducted only recently by Baraffe et 
al. (2003) and Burrows et al. (2003). The calcula- 
tions presented here proceed as described in Chabrier & 
Baraffe (1997) for the non-irradiated case. A grid of ir- 
radiated atmospheric structures is computed for various 
(log Tcft) conditions, for a given orbital distance a, and 
interpolation between these structures yield the unique 
boundary condition with the internal structure, at large 
enough optical depth (t w 100). As shown initially by 
Guillot et al. (1996), irradiation pushes the radiative- 
convective boundary deeper into the interior and, under 
strong irradiation conditions, may force the boundary 
between the atmospheric and internal structure to lie in 
a radiative layer. To insure consistency between the at- 
mosphere and internal thermal structures, we have com- 
puted the Rosseland mean of the atmospheric opacities 
and used this value to calculate the radiative gradient for 
the interior structure. The evolution of the planet obeys 
the usual first and second principles of thermodynamics. 
The energy balance, however, must include the incoming 
stellar radiative flux J-mc = l/2(-R*/a)'^J-"*, where 
and Ri, denote the parent star flux and radius, respec- 
tively. The total emergent flux of the planet J-out now 
reads : 

J' out = ^ inc + o" 

= A T,nc + (1 - ^) T,n, + <JT^S (1) 

where Toff denotes the intrinsic effective temperature of 
the planet, A its Bond albedo, and the first term on 
the r.h.s. of eqn.JQl describes the reflected part of the 
spectrum. Note that the albedo A is not a free parameter 
in our calculations but is calculated consistently with the 
radiative transfer equation. For both transit cases, the 
models predict A < 0.1. The total luminosity thus reads: 

Ltot— Lre fleeted + 47r_Rpf7(Tgg + T^^ff) 

= Lre fleeted + 4:71 R^aTl^^ ~ ( T^dm (2) 



where S denotes the specific entropy of the planet. 
Whereas the total flux J-out or luminosity Ltot is the 
quantity accessible to observation, only the last term 
on the right hand side of eqn.Q concerns the evolution 
of the planet intrinsic luminosity Lint- The first term, 
which defines the Bond albedo, illustrates the fraction of 
the stellar luminosity reflected by the planet atmosphere. 
The second term defines its equilibrium temperature, i.e. 
the temperature the planet would reach in the absence of 
any internal source of energy (Tcff — > 0). Figure |21 com- 
pares the emergent flux calculated for OGLE-TR-56b 
and HD209458b at their present ages. In both cases, the 
infrared part of the spectrum ( ^ 1 ^m) is dominated by 
the re-radiation of the absorbed incident stellar flux. In- 
deed, under the conditions of interest, this contribution 
largely dominates the intrinsic contribution of the planet 
{Teq ~ 2400 K > Teff ~100 K for OGLE-TR-56b, and 
Teq ~ 1700 K > Te// -100 K for HD209458b). For 
IID209458b, the short wavelength part of the spectrum 
(<0.5/im) is mostly due to reflection of stellar light by 
II2 Rayleigh scattering. However, in OGLE-TR-56b the 
high equilibrium temperature, Teq = 2400 K, leads to 
a lower concentration of II2. Consequently, reflection 
is less significant and the majority of the spectrum of 
OGLE-TR-56b is thermal radiation. 

Figure 13 displays the evolution of the radius for dif- 
ferent strongly irradiated planet conditions, with masses 
0.69 Mj„p and 1.5 Mj„p, orbiting a G2 star {T^s = 5900 
K) at different orbital distances, namely a =0.046 AU 
and a =0.023 AU. At a given orbital distance, the less 
massive the planet the larger the effect of irradiation. 
This was expected from Figure Ql showing the larger 
modification of the inner atmosphere profile with larger 
stellar/planet flux ratio. In all cases, however, irradiation 
substantially slows down the contraction of the planet, 
in particular at the early stages of evolution, when de- 
generacy effects in the interior are smaller. The case 
m = 1.5Mj„p and a =0.023 AU (thick solid line) corre- 
sponds to OGLE-TR-56b, for its recently revised mass, 
m ~ 1.45±0.23Mj„p, and radius 1.23±0.15 Rj„p (Torres 
et al. 2003). The age of the system was derived from the 
evolution of the parent star for its observed extinction- 
corrected magnitude and effective temperature (Sasselov 
2003). We derive m ~ 1.05 M©, t ^ A Gyr to reproduce 
the observations, in excellent agreement with Sasselov 
(2003). The striking result is the very good agreement, 
within the error bar, between the present calculations 
and the observations. A prediction of the theory is that 
irradiation from a G2 star yields less than about 15% in- 
crease of the planet radius after 1 Gyr, compared with 
the non-irradiated case, with all radii merging within 
the — 1.1-1.15 Kjup range at this stage. A key test for 
the theory, however, would be the observation of young 
( < 10^ yr) transits, where radii of irradiated planets are 
predicted to vary rapidly with time and to depend sub- 
stantially on the orbital distance. 

The case of IID209458b thus becomes a puzzle. As 
discussed at length in Baraffe et al. (2003), no consis- 
tent evolution calculation can reproduce the observed 
radius. Note that our calculations take into account 
the extension of the atmosphere, which we find to be 
of the order of 0.05 Rj^p (at optical wavelengths) for 
the age of this object. Given the remaining uncertainties 
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in the opacities of these planets, including for example 
non-equilibrium processes or day-night temperature gra- 
dients, we cannot exclude larger extensions. It is very un- 
likely, however, that different opacities modify the inter- 
nal adiabat strongly enough to increase the planet radius 
by 30% and thus solve the aforementioned discrepancy 
between theory and observation for HD209458b. Such a 
modification requires a drastic change of the energy de- 
posited at the top of the internal adiabat, which essen- 
tially determines the radius (Baraffe et al. 2003, Guil- 
lot & Showman 2002). Molecular absorption, however, 
will prevent incident photons from penetrating that deep 
(T ~ 2000 K near the top of the internal adiabat, see 
Figure HJ. As discussed in §4.3 of Baraffe et al. (2003), 
dynamical processes such as tidal dissipation or synchro- 
nization are unlikely to provide the extra source of en- 
ergy, for these timescales are at most of the order of 10* 
yr. An interesting alternative scenario involves the inner 
redistribution of large scale atmospheric kinetic energy 
(winds) generated by the incident irradiation (Guillot 
& Showman 2002, Showman & Guillot 2002). Besides 
arbitrarily modifying the atmosphere-interior boundary 
condition, however, these calculations face a similar issue 
as mentioned above, namely about 1% of this kinetic en- 
ergy must constantly be dissipated at the internal adiabat 
level, i.e. P - 500 bar for HD209458b (see Figure 5 of 
Baraffe et al. 2003) and - 4 x 10^ bar for OGLE-TR-56b 
(see Figure Pi . This corresponds to 18 and 23 pressure 
scale heights for these two planets, respectively. Whether 
such a fraction of the surface flux of kinetic energy can 
be transported downward over so many scale heights and 
converted into thermal energy remains to be uniquely de- 
termined. 

In order to test this scenario in the case of 
OGLE-TR-56b, we have made evolution calculations 
where a fraction of the incident stellar flux is deposited 
deep enough in the planet to modify its internal adi- 
abat (see Baraffe et al. 2003). Figure ^ displays the 
evolution for the two transiting planets when a fraction 
Xi, = Tdep/^inc of the incident flux over the day-side 
surface is deposited in the interior, with Xi, ~ 0.1% to 1 
or 2%. These values correspond to about 40 to 400 times 
the planet's intrinsic energy flux for OGLE-TR-56b and 
up to about 2000 times for IID209458b. As shown in the 
figure, if indeed OGLE-TR-56b is unambiguously con- 
firmed as a transiting planet, values of X^, larger than 
^ 0.5% seem to be excluded, for they yield an interior en- 
tropy and radius that are too large. For HD209458b, on 
the contrary, Xi, = 0.5% is about the minimum fraction 
needed to be consistent with the observational error bar. 
As shown earlier, however, the internal adiabat for the 
IID209458b irradiation conditions lies at fewer pressure 
scale heights from the surface than for OGLE-TR-56b so 
that we expect intuitively a larger fraction of energy to 
be transported downward. Although these conclusions 
must be considered with caution, they suggest that the 
dissipative process invoked by Showman & Guillot may 
indeed take place in hot-jupiters. A tiny ( < 0.5%) frac- 
tion of the surface kinetic flux transported downward and 
converted into thermal energy at the planet's radiative- 
convective boundary is enough to modify significantly 
its internal adiabat and thus its contraction rate. The 



result would be a total ~ 20-30% increase of the radius, 
instead of the aforementioned ~ 15% increase due to bare 
insulation. Another appealing alternative solution to the 
IID209458b dilemma is the suggestion by Bodenheimer, 
Laughlin & Lin (2003) of a second orbiting planet, forcing 
the finite eccentricity (e ~ 0.03) of HD209458b. The con- 
stant tidal heating due to the orbit circularization could 
provide the required mysterious extra source of energy. 

This key issue will be nailed down in the near future 
with further observations of transits, as planned with the 
various ground-based (STARE, VULCAN) and spaced- 
based (COROT, MONS, MOST) missions. 

3. CONCLUSION 

In this letter, we have extended the theory we derived 
recently for the evolution of HD209458b to various con- 
ditions of so-called hot-jupiters, i.e. gaseous planets in 
close enough orbit around their parent star for the incom- 
ing stellar irradiation to affect substantially their struc- 
ture and evolution. As argued previously, a consistent 
treatment between the atmospheric profile and the inte- 
rior profile is mandatory to yield reliable results, both in 
the irradiated and non-irradiated case. The theory, in its 
most simple foundation (ID irradiation, with incoming 
flux deposited uniformly over the day side and no dynam- 
ical redistribution of the incident flux due to day-night 
temperature differences) is found to adequately repro- 
duce the radius of OGLE-TR-56b. These conclusions, of 
course, rely on OGLE-TR-56b being a genuine exoplanet. 
Recent high resolution spectroscopic observations seem 
to support this conclusion (Konacki et al. 2003). If 
OGLE-TR-56b is confirmed to be, unambiguously, a 
planet, then the present results add intriguing support to 
the suggestion that the large radius of IID209458b might 
be due to ongoing tidal interaction with an undetected 
companion. However, we show that a very small fraction 
of the stellar incident flux, transformed into large scale 
kinetic energy partly dissipated at the internal adiabatic 
level, is sufficient to modify substantially a planet's in- 
ternal entropy rate and thus evolution. This fraction is 
unlikely to be larger than a few 0.1%, but certainly de- 
pends on the irradiation conditions. A fraction ^ 0.1 
to 0.5 % would superbly reproduce the observations of 
the two presently detected transits, OGLE-TR-56b and 
IID209458b. The present theory allows us to derive the 
fundamental properties of hot-jupiters (mass, tempera- 
ture, radius) and their observable spectroscopic or pho- 
tometric signatures as a function of age and orbital dis- 
tance. Additional transit detections will also allow bet- 
ter constraints on possible irradiation induced dissipative 
processes, as explored in this letter. Comparisons be- 
tween the present evolution theory and transit or direct 
detections of exoplanets will greatly improve our under- 
standing of these objects and provide strong motivation 
for exploring this exciting new field of astronomy. 

Part of this work was conducted at the Anglo- 
Australian Observatory, under the auspice of the franco- 
australian PICS program. This research was supported 
in part by the LTSA grant NAG 5-3435 to Wichita State 
University. 
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Fig. 1. — Atmospheric profiles (P in cgs units) down to the internal adiabat for non-irradiated (long-dash line) and irradiated planets 
with Tcff = 500 K and T^g = 100 K, logg = 3.5 orbiting a G2 star at different orbital distances: a = 0.046 AU (dash-dot line) and 
a = 0.023 AU (solid line). Symbols indicate various optical depths (defined at A = 1.2/im) and the top of the convective zone. 
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Fig. 2.— The emergent fluxes for T^s = 100 K irradiated models corresponding to OGLE-TR-56b (logg = 3.5) and HD209458b 
(logg = 3.0). The lower two curves that peaJc between 0.3 and 0.4 fim are the reflected fluxes for OGLE-TR-56b (grey) and HD209458b 
(black). The dotted line is the continuation of the thermal flux for HD209458b. All spectra have been reduced to a resolution of 20A. 



Evolution of irradiated planets 



7 




6 7 8 9 10 

log t 

Fig. 3. — Evolution of the radius, in units of Jupiter radii, of irradiated planets with different mass: m=0.69 Mjj,p (thin lines) and m=1.5 
Mjup (thick lines), orbiting a G2 star at different orbital distances: a = 0.023 AU (solid lines) and a = 0.046 AU (dashed lines). The case 
1.5 M 

jup^ ^ — 0.023 corresponds to OGLE-TR-56b. Dash-dotted curves display the evolution in the non- irradiated case, corresponding to 
planets far away from their parent star (a ^ 1 AU). 
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Fig. 4. — Same as Figure|21when various fractions of the incident luminosity Line are included as an extra source of energy in the planet 
interior. The cases are representative of OGLE-TR-56b (upper panel) and HD209458b (0.69 Mj, a=0.046, lower panel). The 0% case (solid 
lines) correspond to the standard irradiated case. Note that Line is ~ 4 times larger in the case of OGLE-TR-56b than for HD209458b. 



